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XMM-NEWTON AND CHANDRA OBSERVATIONS OF THE GALAXY GROUP NGC 5044. 
II. METAL ABUNDANCES AND SUPERNOVA FRACTION 
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ABSTRACT 

Using new XMM and Chandra observations we present an analysis of the metal abundances of the hot 
gas within a radius of 100 kpc of the bright nearby galaxy group NGC 5044. Motivated by the inconsistent 
abundance and temperature determinations obtained by different observers for X-ray groups, we provide 
a detailed investigation of the systematic errors on the derived abundances considering the effects of the 
temperature distribution, calibration, plasma codes, bandwidth, Galactic Ah, and background rate. The 
iron abundance (Zp e ) drops from Z-p e « IZq within R ss 50 kpc to Zp c « 0.4 solar near R = 100 kpc. This 
radial decline in Zp c is highly significant: Zp c = 1.09±0.04Z© (statistical) ±O.O5Z +0.18Z Q (systematic) 
within R = 48 kpc (5') compared to Z Fc = 0.44 ± O.O2Z (statistical) ± 0.10Z Q + 0.13Z Q (systematic) 
over R = 48 — 96 kpc (5' — 10') . There is no evidence that the radial profile of Zy Q flattens at large radius. 
The data rule out with high confidence a very sub-solar value for Z-p e within R — 48 kpc confirming 
that previous claims of very sub-solar central Zp c values in NGC 5044 were primarily the result of 
the Fe Bias: i.e., the incorrect assumption of spatially isothermal and single-phase gas when in fact 
temperature variations exist. Next to iron the data provide the best constraints on the silicon and sulfur 
abundances. Within R = 48 kpc we obtain Z si / Z Fc = 0.83 ± 0.02 (statistical) ± 0.02 + 0.07 (systematic) 
and Zs/Zp c = 0.54 ± 0.02 (statistical) ± 0.01 + 0.01 (systematic) in solar units. These ratios (1) are 
consistent with their values at larger radii, (2) strongly favor convective deflagration models over delayed- 
detonation models of SNe la and (3) imply that SNe la have contributed 70%-80% of the iron mass within 
a 100 kpc radius of NGC 5044. This SNe la fraction is also similar to that inferred for the Sun and 
therefore suggests a stellar initial mass function similar to that of the Milky Way. We mention that at 
the very center (R « 2 kpc) the XMM and Chandra CCDs and the XMM RGS show that the Fe and Si 
abundances drop to « 50% of their values at immediately larger radius analogously to that seen in some 
galaxy clusters observed with Chandra. We find the magnitude of this dip to be sensitive to assumptions 
in the spectral model, but if real it is difficult to reconcile with the expectation that metal enrichment 
from the stars in the central galaxy should result in a centrally peaked metal abundance profile in the 
hot gas. 

Subject headings: X-rays: galaxies: clusters - galaxies: halos - galaxies: formation - cooling flows - 
galaxies: individual: NGC 5044 
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1. INTRODUCTION 

There is presently a controversy associated with the iron 
abundances of groups (and the most X-ray luminous el- 
liptical galaxies) deduced from X-ray observations. While 
there seems to be general agreement of sub-solar iron abun- 
dances outside the central regions (r w 50 — 100 kpc) of 
groups (e.g., Finoguenov & Ponman 1999; Buote 2000a), 
different investigators have often obtained (for the same 
groups) different results for the central regions (r < 
50 kpc) where the metal enrichment from a central galaxy 
should be most pronounced. Most previous ROSAT and 
ASCA studies have found very sub-solar values of Z-p e in 
the central regions of groups (for reviews see, Buote 2000b; 
Mulchaey 2000). Since these low values of Zp c are gener- 
ally lower than the stellar iron abundances (e.g., Trager 
et al. 2000), they imply that Type la supernovae (SNe la) 
cannot have contributed significantly to the enrichment of 
the hot gas. This implies that there is a lower binary star 
fraction and SNe la rate in the group galaxies so that most 
of the iron derives from SNe II with a "top heavy" stel- 
lar initial mass function (IMF) (e.g., Renzini et al. 1993; 



Renzini 1997; Arimoto et al. 1997). Consequently, various 
authors have questioned the reliability of X-ray determi- 
nations of Zp c and have suggested that the low Zp e values 
are caused by errors associated with the Fe L lines in X-ray 
plasma codes (e.g., Arimoto et al. 1997; Renzini 2000). 

However, in a series of papers (Buote & Fabian 1998; 
Buote 1999, 2000b, a) we found that indeed the iron abun- 
dances in the central regions of groups were measured in- 
correctly, but not because of errors in the plasma codes. 
Instead, we attributed the very sub-solar Zp c values to 
an "Fe Bias" arising from forcing a single-temperature 
model to fit a spectrum consisting of multiple tempera- 
ture components with temperatures near 1 keV (see es- 
pecially, Buote 2000b, a). The multiple temperature com- 
ponents may arise either from a radially varying single- 
phase gas or represent real multiphase structure in the 
hot gas. We found near-solar values for Zp c within the 
central 50-100 kpc of groups, which is larger than the typ- 
ical stellar value for Zp e within R e /2 in elliptical galaxies 
(Trager et al. 2000), implying that a significant number 
of SNe la have enriched the hot gas, in better agreement 
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with a Galactic IMF. 

Even stronger constraints on the SNe la fraction and 
the IMF are placed by the ratios of the abundances of a 
elements to iron (e.g., Gibson et al. 1997; Renzini 1997; 
Brighenti & Mathews 1999). Previous X-ray observations 
did not place strong constraints on the a abundances in 
groups and were usually consistent with solar values (e.g., 
Mulchacy 2000, and references therein). 

Recently, using a new XMM observation of the luminous 
X-ray group NGC 1399 we find Z Fc /Z Q w 1.5 - 2 within 
r ~ 20 kpc and Zsi/^Fe ~ 0.8 solar over the entire 50 kpc 
radius studied. These results imply a SNe la fraction of 
« 80% which is similar to that inferred for the Sun and 
therefore suggests a stellar initial mass function similar to 
the Milky Way as advocated by Renzini and others (e.g., 
Renzini et al. 1993; Renzini 1997; Wyse 1997). 

NGC 5044 is brighter and more luminous in X-rays than 
NGC 1399 but is slightly lower in temperature. In Buote 
et al. (2003, hereafter Paper 1) we showed using XMM 
and Chandra data that within r « 30 kpc the hot gas is 
not isothermal, nor is it consistent with a radially vary- 
ing single-phase medium. Instead a limited multiphase 
medium is implied where the temperature varies from ~ Th 
to - T h /2 (T h w 1.3 keV), but no lower. 

In this paper we measure the metal abundances of the 
hot gas in NGC 5044 using the XMM and Chandra data 
and provide a detailed investigation of the systematic er- 
rors in the derived abundances considering the effects of 
the temperature distribution, calibration, plasma codes, 
bandwidth, Galactic TVh, and background rate. The impli- 
cations of these measurements for the supernova fraction 
and IMF are then discussed. 

The paper is organized as follows. In §2 we present the 
iron abundance as a function of radius for different spectral 
models. We present the profiles of silicon in §3 and other 
abundances in §4. A comprehensive discussion of system- 
atic errors in the abundance measurements is given in §6. 
We give in §7 our most precise constraints for the emission- 
weighted average abundances in regions of 0-50 kpc and 
50-100 kpc with a full and explicit accounting of the rele- 
vant systematic errors. Finally, in §8 we present our con- 
clusions. We assume a distance to NGC 5044 of 33 Mpc 
using the results of Tonry et al. (2001) for Ho = 70 km 
s" 1 Mpc" 1 (note: 1" = 0.160 kpc). 

2. IRON ABUNDANCE 

2.1. Preliminaries 

To obtain the three-dimensional properties of the X-ray 
emitting gas we perform a spectral deprojection analysis 
assuming spherical symmetry using the (non-parametric) 
"onion-peeling" technique as discussed in Paper 1. We re- 
fer to deprojected models as "3D" while traditional model 
fits to the data on the sky are referred to as "2D" (i.e., 
without deprojection). However, with respect to 2D mod- 
els, deprojection always inflates the errors between points 
which is related to the error associated with the deriva- 
tive of the emissivity in an Abel inversion. As discussed 
in Paper 1, we perform a rcgularization procedure on the 
oxygen and neon abundances to limit their radial fluctua- 
tions in 3D models. No rcgularization is applied to any 2D 
model. Because the abundances obtained from 2D models 
have smaller statistical errors and do not require any reg- 



ularization we shall generally present results for both 2D 
and 3D models in this paper. Statistical errors on the pa- 
rameters are computed using the Monte Carlo procedure 
discussed in §4.1 of Paper 1. 

We measure the Fe abundance as a function of radius 
using a suite of different models for the temperature distri- 
bution as described in Paper 1. Single-temperature (IT) 
and two-temperature (2T) models are used as our baseline 
models for comparison to previous studies. We also exam- 
ined a set of models that emit over a continuous range 
of temperatures; i.e., models having a continuous differen- 
tial emission measure (DEM): cooling flow, gaussian DEM, 
and power-law DEM (PLDEM) . In every model the follow- 
ing abundances are free parameters: Fe, O, Ne, Mg, Si, 
and S - the abundances for all other elements are tied to 
Fe in their solar ratios. Unless stated otherwise, for every 
multitemperature model (e.g., 2T) the abundances of one 
temperature component are tied to those in the other tem- 
perature component (s). Even though the 2T and PLDEM 
models provide superior fits within the central ~ 30 kpc of 
NGC 5044, in this paper we present results for all models 
for the purpose of showing the dependence of the inferred 
Fe abundance on the temperature structure of the hot gas. 

The reference solar Fe abundance has been a source 
of much confusion in the literature. There is now good 
agreement between values of the Fe abundance obtained 
from measurements in the solar photosphere and from 
solar-system meteorites (e.g., McWilliam 1997; Grevesse 
& Sauval 1998). Therefore, we take the solar abundances 
in XSPEC (vll.2.0af) to be those given by the Grevesse & 
Sauval (1998) table which use the correct Fe value, Fe/H = 
3.2 x 10 by number. Unfortunately, most previous and 
many current X-ray studies of Fe abundances use the in- 
correct "old photospheric" value of Fe/H = 4.7 x 10~ 5 
present in the Anders & Grevesse (1989) table in xspec 
which is still the default option in xspec. Consequently, 
investigators who use the old photospheric value for Fe/H 
obtain values for the Fe abundance that are approximately 
a factor 1.4 too small. In comparing with our results, we 
shall transform all abundances from previous studies to 
those of Grevesse & Sauval (1998) unless stated otherwise. 

2.2. Results 

The Fe abundance (Z Fc ) as a function of radius obtained 
from IT models fitted simultaneously to the XMM and 
Chandra data is displayed in Figure 1 (Left panel); note 
that the Chandra data only apply to the inner three ra- 
dial bins. In Table 1, for convenience, we reproduce the 
temperature results for the IT and 2T models from Pa- 
per 1. For r < 50 kpc the iron abundance is approx- 
imately constant such that Z Fc ~ 0.65Z Q for the IT 
(2D) model and Z Fc « 0.75Z Q for the IT (3D) model. 
At larger radii Z Fe decreases with increasing radius with 
the lowest value occurring in the bounding annulus/shell: 
Z Fc = 0.28 ± 0.02Z Q for r = 73 - 97 kpc. The ~ 15% 
larger values of Z Fe obtained for the deprojected models 
near 40-50 kpc arc primarily the result of the projection 
effect which tends to smear out a radially varying func- 
tion, while for smaller radii they are primarily the result 
of the Fe Bias. In the latter case the deprojection removes 
the projected temperature components from exterior shells 
thus allowing the IT model to be a better - though still 



3 



10 



50 100 



50 100 



s 



1 1 1 1 
■ 1 T 


i i 1 1 i i i | i i i i | 

- 


- 


3D 




" □ ""^-"••'"n'S - ( D / N ' 


. i . 


2D : 

■ ■ I ft*' 



0.5 1 

r (arcmin) 




0.4 



2T 



0.5 1 

r (arcmin) 



Fig. 1. — Radial profiles (units — bottom: arcmimites, top: kpc) of the Fe abundance and associated lc errors for (Left panel) IT and 
(Right panel) 2T models fitted simultaneously to the XMM and Chandra data. Note that the Chandra data apply only to the inner three 
radial bins. In each case "3D" refers to results obtained from a spectral deprojection analysis of the XMM and Chandra data. 




not good - representation of the spectrum within a given 
shell at smaller radii. 

Previous ASCA studies of NGC 5044 that fitted IT 
models to the w 0.5 — 10 keV data in the central re- 
gions obtained much smaller abundances. For example, 
Fukazawa et al. (1996) obtained Z Fe = 0.46 ± 0.03Z Q for 
R < 20 kpc, and Finoguenov & Ponman (1999) obtained 
Z Fc = O.65±O.O7Z0 for r = 30-40 kpc using IT models of 
the ASCA data of NGC 5044. As we discussed in Buote 
(1999), IT models are poor fits to the ASCA data ex- 
tracted from large apertures because they contain a distri- 
bution of temperatures implied by the radial temperature 
gradients observed with ROSAT (David et al. 1994; Buote 
1999, 2000a). These temperature gradients suggested by 
ROSAT are confirmed and mapped in detail with Chandra 
and XMM in Paper 1. 

The Fe abundances for the 2T models fitted simultane- 



ously to the XMM and Chandra data are also displayed 
in Figure 1 (Right panel). For all radial bins interior to 
a radius of « 40 kpc, excluding the central bin, the val- 
ues of Z Fc obtained from the 2T models exceed by 40%- 
100% those obtained from the IT models. This system- 
atic increase arises from the Fe Bias, and the differences 
between the IT and 2T models are highly significant; e.g., 
in shell #2 we obtain Z Fc = 0.78 ± O.O3Z for IT (3D) 
Z Fc = 1.24±0.10Z Q for 2T (3D), and Z Fc = 1.18±O.12Z 
for PLDEM (3D) models. Recall from Paper 1 that the 
2T and PLDEM models provide the best fits to the data, 
and the fit residuals observed for the IT model near 1 keV 
are fully consistent with those produced by the Fe Bias. 

In Figure 2 we show the Fe abundances obtained for 
the PLDEM model. It is notable that the PLDEM model, 
which fits the data almost as well as the 2T model, gives 
Z Fe values consistent with the 2T model everywhere within 
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Table 1 

Temperatures for IT and 2T Models 









IT (2D) 


IT (3D) 


2T (2D) 


2T (3D) 






Rout 


T 


T 


T c 


T h 


T c 


T h 


Bin 


(arcmin) 


(arcmin) 


(kcV) 


(kcV) 


(keV) 


(keV) 


(keV) 


(keV) 


1 


0. 


0.5 


0.754 ± 0.005 


0.698 ±0.013 


0.74 ±0.01 


1.66 ±0.09 


0.69 ±0.01 


3.7(> 1.4) 


2 


0.5 


1.5 


0.807 ± 0.002 


0.788 ± 0.003 


0.787 ±0.003 


1.47 ±0.04 


0.772 ± 0.004 


1.62 ±0.07 


3 


1.5 


2.5 


0.985 ± 0.003 


0.904 ± 0.006 


0.83 ±0.01 


1.35 ±0.02 


0.82 ±0.01 


1.36 ±0.06 


4 


2.5 


3.5 


1.172 ±0.005 


1.111 ±0.018 


0.87 ±0.01 


1.36 ±0.02 


0.86 ±0.01 


1.37 ±0.06 


5 


3.5 


4.5 


1.236 ±0.010 


1.251 ±0.019 


0.92 ±0.01 


1.31 ±0.02 


1.05 ±0.08 


1.47 ±0.13 


6 


4.5 


5.7 


1.229 ±0.007 


1.244 ±0.017 


0.85 ±0.10 


1.35 ±0.03 


0.86 ±0.14 


1.37 ±0.04 


7 


5.7 


7.6 


1.205 ±0.011 


1.268 ±0.025 


0.79 ±0.14 


1.29 ±0.05 


0.50 ±0.23 


1.28 ±0.05 


8 


7.6 


10.1 


1.147 ±0.013 


1.147 ±0.015 


0.90 ±0.07 


1.30 ±0.04 


0.90 ±0.10 


1.30 ±0.05 



Note. — Temperature results obtained in Paper 1 for IT and 2T models. These models have been fitted simultaneously to the XMM and 
Chandra data in bins 1-3 and only to the XMM data for bins 4-8. "Bin" refers either to a circular annulus (2D) or spherical shell (3D). For 
the case where a lower limit is given it represents the lowest value obtained from 100 error simulations; i.e., it is essentially a 99% confidence 
lower limit. 



the lcr — la errors. The overall similarity of Z Fc values for 
different multitemperature models (including cooling flow 
and Gaussian DEM - not shown) indicates that the value 
of Z Fo can be fairly reliably inferred with any of the mul- 
titemperature models; i.e., it is principally the IT model 
that suffers from the Fe Bias and gives large underesti- 
mates of Z Fc . 

Focusing on the multitemperature models since they 
provide the best spectral fits in the central regions, we con- 
clude that the iron abundance displays a strong gradient in 
NGC 5044. At the largest radius examined (R = 97 kpc) 
we have Z Fe w 0.35Z Q . The (deprojected) iron abundance 
increases to values between f-f .5 solar within r « 50 kpc 
for all multitemperature models and then dips within the 
central bin to a value of ~ O.8Z . 

For every model we have discussed so far (except the 1 T 
(2D) model) the iron abundance dips in the central bin. 
In §8 we discuss a multi-temperature model that does not 
dip in the center and the implications of such a dip for 
enrichment models of the hot gas. 

3. SILICON ABUNDANCE 

Next to the broad feature of Fe L lines near 1 keV, the 
most notable emission lines in the EPIC spectra of NGC 
5044 (see Paper 1) are the Ka lines of silicon; i.e., Si XIII 
Hea (1.85 keV), Si XIV Lya (2.0 keV). It follows that the 
silicon abundance (Zsi) is also best constrained next to 
Z Fe . Although it might be expected that Z$i does not suf- 
fer the same model-dependences as Z Fe because the silicon 
lines are fairly isolated and well-separated from the Fe L 
lines, the results of the spectral fits show otherwise. 

In Figure 3 we display Zgi as a function of radius for 
IT and 2T models. Analogously to Z Fe we find that Z$i 
obtained from deprojected fits (i.e., 3D) of both IT and 
2T models arc systematically larger than those obtained 
from 2D fits by 10%-20%. For IT models, Z si w Q.7Z Q 
for r < 30 kpc and falls to Z$i ~ 0.2Z Q in the outermost 
radial bin. Similar values are obtained at large radius for 
2T models. But for r < 30 kpc, 2T (3D) model yields 
Z si w 1.1Z for r < 30 kpc which dips to Z si w 0.7Z Q 
in the central bin in good agreement with the behavior 
observed for Z Fc . 

The other multitemperature models (cooling flow, Gaus- 



sian DEM, power law DEM) give values for Zgi relative to 
the IT and 2T models that are entirely analogous to that 
described for Z Fe in the previous section; i.e., they mostly 
give values consistent with the 2T model within the 1-2 a 
errors. Given that the 2T and PLDEM models provide the 
best fits within r < 30 kpc (Paper 1), and they yield fully 
consistent values for Z Fe and Z$i, the near-solar values ob- 
tained for each element should be considered the favored 
values. 

Although the silicon abundance obtained from IT mod- 
els within the central regions is biased low because of mul- 
tiple temperature components having values near 1 keV 
within each radial bin (i.e., "Silicon Bias", Buote 2000b), 
the ratio of silicon-to-fe (Zgi/Z Fe ) is affected much less. 
In Figure 3 (Right panel) we plot Zsi/Z Fe as function of 
radius for IT (2D) and 2T (2D) models; the 3D models are 
everywhere consistent with their 2D counterparts within 
the lcr — 2cr errors. The Gaussian- weighted mean of all 
radial bins is Zgi/ Z Fc = 0.96 ± 0.02 (in solar units) for IT 
(2D) and Z si /Z Fc = 0.83±0.01 (in solar units) for 2T (2D); 
the gas-mass- weighted values are Zsi/Z Fc = 0.74 ± 0.05 
and Z si /Z Fc = 0.74 ± 0.05 respectively for the IT and 
2T models. As expected, the largest differences between 
IT and 2T models are observed within the central four 
bins (i.e., r < 30 kpc) where Z si /Z Fc = 1.00 ± 0.02 and 
Z S i/Z Fc = 0.85 ±0.01 respectively for the IT and 2T mod- 
els for Gaussian means and Zgi/ Z Fc = 0.86 ± 0.03 and 
Zsi/Z Fe = 0.80 ±0.02 respectively for corresponding mass- 
weighted means; for comparison, in the outer bins (5-8) we 
have Z si /Z Fc = 0.74±0.04 (IT) and Z si /Z Fc = 0.76±0.03 
(2T) for Gaussian means and Z Si /Z Fc = 0.72 ± 0.06 (IT) 
and Zsi/ Z Fe = 0.73 ± 0.06 (2T) for mass- weighted means. 

We note that the marginal evidence that Zsi/Z Fc de- 
creases with radius in NGC 5044 (i.e., for the 2T model 
the Gaussian mean value of Zsi/ Z Fe for bins 5-8 is lower 
than that obtained from bins 1-4 by 1.9cr) should be con- 
sidered quite tentative because this relatively small dif- 
ference could be an artifact of the simple spectral models 
used. (However, all models we have investigated give sim- 
ilar results.) If real, the slightly larger a/Fe abundance 
ratios near the center might be the result of stellar mass 
loss which would be expected to have a lower SNe la frac- 
tion. 



5 



10 



50 100 



: IT 


1 , , 1 1 


, , iii' 


- 


3D 


- 




2D 


JL >\ ' " □ ) 

<3>/ a \; D .\ - : 


1 I 


, , , 


<•>: 

■ i i 1 1 , , 




r (arcmin) 



50 100 



50 100 



2T 



_i 




r (arcmin) 



r (arcmin) 



Fig. 3. — Radial profiles (units — bottom: arcminutes, top: kpc) of the Si abundance and associated la errors for (Left panel) IT and 
(Middle panel) 2T models. In each case "3D" refers to results obtained from a spectral deprojection analysis. (Right panel) Zsi/Z-pe ratios 
for the IT (2D) and 2T (2D) models. 



4. OTHER ABUNDANCES 

In Figure 4 we show the radial abundance profiles for O, 
Mg, and S for the IT (2D) and 2T (2D) models; the de- 
projected values (not shown) are considerably noisier and 
are in most cases consistent with the 2D values within the 
la errors. Although there is considerable scatter between 
radii, like Si the values of the O and Mg abundances ob- 
tained with the 2T model are systematically larger than 
those obtained with the IT model within r ~ 50 kpc. The 
values of S are generally consistent between the IT and 2T 
models. The relative insensitivity of S to the temperature 
model is likely attributed to the fact that the key S emis- 
sion lines near 2.4 keV are farther away from the ~ 1 keV 
Fe L lines than either O, Mg, or Si. 

Although the value of the Mg abundances differs sig- 
nificantly between IT and 2T models, the Mg/Fe ratio 
- like the Si/Fe ratio - is quite similar for both models; 
e.g., in bin #2, Zm s / Z Fe = 0.94 ± 0.05 in solar units, for 
IT (3D) while, Z Mg /Z Fc = 0.88 ± 0.05 in solar units, for 
2T (3D). But the 6/Fe and S/Fe ratios do show signifi- 
cant differences between the models, although the O/Fe 
ratio is very sub-solar in both cases; e.g., in bin #2, 
Z /Z Fc = 0.19 ± 0.04 and Z s /Z Fc = 0.96 ± 0.08 in so- 
lar units, for IT (3D) while, Z /Z Fc = 0.30 ± 0.05 and 
Z s /Z Fc = 0.66 ± 0.05 in solar units, for 2T (3D). We note 
that the oxygen abundance in bin #7 is significantly larger 
(« 4<r) than its values in adjacent bins. We do not believe 
the fluctuation in the oxygen abundance in bin #7 is phys- 
ical, but we have not yet isolated the cause. 

For the 2T (3D) model we obtain Gaussian weighted 
mean values of all radial bins of Zs/Z Fc = 0.62 ± 0.04, 
Z /Z Fc = 0.39 ± 0.03, and Z Mg /Z Fc = 0.88 ± 0.04 in 
solar units; the corresponding gas-mass- weighted means 
are Z s /Z Fc = 0.46 ± 0.10, Z /Z Fc = 0.56 ± 0.09, and 
Zwig/ Z Fc = 0.66 ±0.17 in solar units. Ignoring the anoma- 
lous bin #7 we obtain Zo/Z Fc = 0.37 ± 0.03 (Gaussian 
averaged) and Zo/Z Fc = 0.37 ± 0.04 (mass averaged) in- 
cluding only bins 1-6. 

Finally, we mention that the Ne abundance (not shown) 
is the most poorly constrained abundance that we have in- 
vestigated. Undoubtedly the poor constraints are at least 



partially attributed to the location of the strongest Ne 
lines near 1 keV which are therefore completely blended 
with the Fe L lines in the CCD spectra. Generally, we find 
very sub-solar values for the Ne abundance except within 
r ~ 10 kpc where near solar values are suggested. More 
precise constraints are provided for the emission-weighted 
average Ne abundance in §7. 

5. NON-AZIMUTHALLY SYMMETRIC ANALYSIS 

Following our discussion in §5 of Paper 1 we searched 
for azimuthal variations in the iron abundance by fitting 
IT and 2T models to spectra extracted from circular aper- 
tures arranged in a square array on the EPIC MOS images. 
Consistent with the results for the temperatures discussed 
in Paper 1, we find that overall the abundances obtained 
from this analysis are consistent with the spherically sym- 
metric analysis. In particular, we find no evidence for 
azimuthal abundance variations associated with the small 
azimuthal temperature variations between radii of 2' — 3' 
mentioned in Paper 1. We also observe the same increase 
in the iron abundance for apertures within r ~ 30 kpc 
for 2T over IT models that is obtained in the spherically 
symmetric (i.e., azimuthally averaged) analysis. 

6. SYSTEMATIC ERRORS 

This section contains a detailed investigation of system- 
atic errors on the abundance measurements. Those readers 
who are not interested in these technical details can safely 
skip ahead to §7. 

6.1. Calibration 

We have examined possible systematic errors in the mea- 
surements of the metal abundances arising from calibra- 
tion differences between the XMM MOS, XMM pn, and 
Chandra ACIS-S3 CCDs. In Table 2 we list the values 
of Z Fc obtained from IT (2D) and 2T (2D) models fitted 
separately to the XMM and Chandra data; i.e., the MOS 
and pn data were fitted simultaneously while the Chandra 
data were fitted alone. We focus on 2D models so that 
the fits for a specific annulus are independent of results 
obtained from fits to adjacent regions at larger radii. 
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Fig. 4. — Radial profiles (units - bottom: arcminutes, top: kpc) of the O, Mg, and Si abundances and their associated la errors for the 
(Left panel) IT (2D) and (Right panel) 2T (2D) models. 



Table 2 

Comparison of Fe Abundances from Chandra and XMM 







IT 








2T 










XMM 


Chandra 


% 


XMM 


Chandra 


% 


% 


Annulus 


(free) 


(free) 


(free) 


(free) 


(fix) 


(free) 


(fix) 


(free) 


(fix) 


1 


0.64 ± 0.05 


0.61 ±0.05 


5±1 


1.17 ±0.13 


1.04 ±0.13 


0.83 ± 0.09 


0.90 ±0.09 


29 ± 14 


13 ± 16 


2 


0.58 ±0.02 


0.68 ± 0.03 


-17 ±6 


1.27 ±0.07 


1.20 ±0.07 


1.11 ±0.07 


1.20 ±0.07 


13 ±8 


0±8 


3 


0.65 ± 0.03 


0.67 ±0.03 


-3 ±6 


1.30 ±0.07 


1.24 ±0.07 


0.99 ± 0.08 


1.13 ±0.08 


24 ±9 


9±9 



Note. — All models are 2D to allow an independent comparison of the data sets in the central regions. The iron abundance is expressed 
in solar units. "%" is the percent difference between the XMM and Chandra abundances, "fix" means that the relative normalizations of the 
temperature components for the 2T models are required to have their best-fitting values obtained from the joint XMM and Chandra fits to 
guarantee a consistent comparison. 
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Fig. 5. — Radial profiles (units — bottom: arcminutes, top: kpc) of the Fe abundance and associated ltr errors for IT (blue circles and 
solid diamonds) and 2T (red boxes and dashed diamonds) models. Only the Chandra data are used and each model is 2D; i.e., not obtained 
from the spectral deprojection analysis. We also show for each model (green crosses and dashed diamonds) the results obtained in the wider 
apertures used in the joint XMM and Chandra fits; the results for these fits are given in Table 2. 



Using the IT model the XMM and Chandra data in 



annuli 1 and 3 give values of Zp c that agree within 5%. 
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In annulus 2 the value obtained from XMM is 17% ± 6% 
lower than that obtained from Chandra. The significance 
of this discrepancy is not very high. However, it is note- 
worthy that in annulus 2 we find the largest improvement 
for a 2T model over a IT model (Paper 1). Given the very 
poor fit of the IT (2D) model within annulus 2 the modest 
discrepancy of 17% ± 6% between detectors could simply 
arise from the different energy-dependent sensitivities of 
the detectors. 

The effects of different sensitivities are illustrated by the 
2T (2D) fits also listed in Table 2. If the relative normal- 
izations of the two temperature components are allowed to 
be different for XMM and Chandra (i.e., denoted by "free" 
in Table 2), then we find that the values of Zp e obtained 
from XMM exceed those obtained by Chandra by 15%- 
30%. However, if we fix the relative normalizations of the 
cooler and hotter components to their best-fitting values 
obtained from the joint XMM-Chandra fits (i.e., denoted 
by "fix" in Table 2), then we find no significant differences 
between the XMM and Chandra fits; i.e., the percentage 
differences are all less than the (relatively large) la errors. 
Since the difference of the quality of the fits (not shown) 
is quite negligible between the fixed and free cases, the 
different values of Zp e obtained in the two cases are prob- 
ably attributed to the different energy-dependent sensitiv- 
ities of the XMM and Chandra detectors. We conclude 
that errors in the relative calibration between XMM and 
Chandra contribute errors in the measured value of Zp e of 
less than 20%. 

The agreement of Zp e obtained from the Chandra and 
XMM CCDs also indicates that the larger PSF of XMM 
does not significantly affect our results; i.e., we have chosen 
bins that are sufficiently wide for the XMM PSF - consis- 
tent with our findings regarding the temperature profile in 
Paper 1. 

(We mention that we have performed an identical study 
to assess differences in Z-p e measured between the XMM 
MOS and pn CCDs. For annuli 5-8 we find that the values 
of Zp c obtained from separate fitting of the MOS and pn 
data generally agree within their 1-2 a errors. However, 
for annuli 1-4 the values of Zp c obtained from the pn data 
always exceed the values obtained from the MOS by as 
little as 10% to as much as 50%, with the MOS values 
usually agreeing better with the Chandra data in their re- 
gions of overlap. As a result, we explored using only the 
MOS and Chandra data in our analysis but found that 
the errors in the deprojected Zp e abundances were suffi- 
ciently large that we had to regularize the iron abundance 
to obtain a radial profile similar to the 2D result. After 
this regularization we noticed that the results were quite 
consistent with the original un-regularized Zp e values ob- 
tained from the simultaneous MOS+pn fits. Because of 
this agreement, and the fact that the MOS+pn fits do not 
require Zp e to be regularized, the good agreement of the 
MOS+pn fits with the Chandra fits discussed above, and 
the excellent agreement of the MOS+pn and ASCA data 
discussed in §7, we decided to use all of the data sets in 
our default analysis.) 

As in Paper 1 we also performed fits to the Chandra 
data within annuli that are half the size of those used for 
the XMM and joint XMM - Chandra fits. In Figure 5 we 
show Z Fc obtained from IT (2D) and 2T (2D) fits to the 



Chandra data in annuli of 30" width between 0.5' — 2.5' 
and 15" width between 0" — 30" (see also Table 5 and 
Figure 6 in Paper 1). Also shown in Figure 5 are the 
results using the larger 1 '-width annuli for the results ob- 
tained from the joint XMM - Chandra fits. The values 
for Zp e obtained in the larger bins are mostly consistent 
with an average value of Zp e obtained in the thinner bins. 
Moreover, if the hot gas is actually single-phase, we would 
have expected in the smaller bins to have a smaller ra- 
dial temperature variation, and thus a weaker "Fe Bias", 
and therefore a smaller underestimate of Zp e using the IT 
model. Instead, we see good agreement between the Zp e 
values obtained using the thinner and wider bins. (Note: 
consistent results are obtained with 3D models using the 
XMM data to deproject shells exterior to r = 2.5'.) This 
agreement is consistent with the models for limited multi- 
phase gas (i.e., particularly the 2T and PLDEM models.) 

In Figure 5 one also notices that the central dip in Zp c 
is apparently a smooth transition beginning within the 
Chandra annulus R — 15" — 30" and finishing within the 
central bin. It is within the central bin (R = 15") where 
the value of Zp e obtained from the IT and 2T models (2D) 
agree best. 

As discussed in §6.1.2 of Paper 1 the value of Zp e ob- 
tained from the XMM RGS data of NGC 5044 by Tamura 
et al. (2003) is in excellent agreement with the value we ob- 
tain from the simultaneous fits to the XMM and Chandra 
CCDs in the central radial bin. We emphasize that this 
excellent agreement occurs when the same models (with 
the same free parameters) are fitted over the same en- 
ergy range for both the RGS and the CCDs. The good 
agreement provides a further check on possible calibration 
differences between detectors. 

We conclude that both the XMM and Chandra CCD 
data require Zp c to be at least solar within r ~ 30 kpc 
except for within the very central region (r = 2.4 kpc). 

6.2. Plasma Codes 

For a single isolated emission line of a particular element 
in a coronal plasma the abundance of that element may 
be directly measured by calculating the ratio of the flux 
within the line to the flux of the local continuum. In such 
an idealized case, any error in the theoretical calculation 
of the line flux within the plasma code will translate di- 
rectly to an error in the measured elemental abundance. 
Consider, however, the case of the Fe L shell lines for a 
system like NGC 5044. 

In a coronal plasma with T ~ 1 keV there is a forest 
of Fe L lines near 1 keV, with the most important lines 
spanning the approximate range 0.7-1.2 keV. Consider also 
these lines when observed at the moderate spectral resolu- 
tion of the XMM and Chandra CCDs (AE w 50 - 100 eV 
near 1 keV). In these cases the Fe abundance is obtained 
by fitting data over a broad energy range, so that all of the 
strong Fe L lines over 0.7-1.2 keV play a key role as does 
the continuum determined from regions outside the Fe L 
region. The net effect is that even rather large errors in a 
small number of lines do not necessarily translate directly 
to large errors in the inferred Fe abundance. 

We illustrate this effect in Figure 6 using the APEC and 
mekal plasma codes. These codes have many differences 
in how they model the atomic physics and in their line li- 
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Fig. 6. — Comparison of the APEC and MEKAL plasma codes in the energy region, 0.7-1.2 keV, where the Fe L shell emission lines dominate 
for a ~ 1 keV plasma. We show the APEC/mekal ratio for a T = 1 keV, solar metallicity plasma (Left paneJ) at high resolution and (Right 
panel) at moderate resolution, the latter corresponding approximately to the response matrix binning of the MOS data of NGC 5044 we used 
in our spectral fitting. 



braries (especially for Fe L) thus allowing for an interesting 
test of the robustness of the iron abundance determination. 
We consider a fiducial spectral model with T = 1 keV, 
Z = IZq, and equal normalizations for the APEC and 
mekal codes. The ratio of the fluxes in the two mod- 
els are plotted as a function of energy over the range 0.7- 
1.2 keV for two different resolutions. In the left panel we 
plot this ratio at high-resolution binning (AE — 0.5 eV) 
giving 1000 energy bins over the energy range shown. Al- 
though most of the points cluster near 1, there are many 
points having large errors; i.e., ratios larger than 2 or less 
than 1/2. In the right panel we plot the same ratio at the 
binning for the moderate resolution CCDs (AE = 15 eV) 
used in our analysis of NGC 5044; specifically, these en- 
ergy bin sizes correspond to those of the re-binned MOS1 
data in the central radial bin for NGC 5044. In this case 
the ratios are all within w 20%, and thus we should expect 
that the iron abundances should not differ by much more 
than this amount for each model. 

Indeed, we find that both IT and 2T models (each 2D) 
computed using the APEC code differ by no more than 10%- 
20% from those computed using the MEKAL code. Inter- 
estingly, the small differences between APEC and mekal 
do appear to be a real systematic effect despite the fact 
that these differences are consistent within the la errors. 
The sign of the systematic difference is opposite for IT 
and 2T models: for IT models, Z Fc apoc > Z Fc mokal , while 
for 2T models, Z Fo apoc < Z Fe mekal . 

We conclude that the measured iron abundances are 
accurate to within 10%-20% considering reasonable re- 
maining errors in the plasma codes such as those asso- 
ciated with the Fe XVII lines near 0.7 keV discussed by 
Bchar et al. (2001). Of course, for high-resolution stud- 
ies using the gratings on Chandra and XMM more care 
must be taken when analyzing the properties of individ- 
ual lines such as done by Behar et al. (2001). But even 
high-resolution studies of NGC 5044 using the RGS obtain 
fully consistent values for Zp e from broad-band spectral fit- 



ting with the values obtained from the XMM and Chandra 
CCDs as discussed at the end of §6.1. 

6.3. Bandwidth 

In section 5.1 of Buote (2000a) we discussed the sensi- 
tivity of the Fe Bias to the lower energy limit used in the 
spectral fitting of a T ~ 1 keV coronal plasma. Here we 
summarize this discussion and refer the interested reader 
to the previous paper for details. Consider (1) a IT model 
(T = 1 keV) of a coronal plasma with solar metallicity, 
and (2) a multitcmperature coronal model with temper- 
ature components distributed across T « 0.7 — 1.5 keV 
also with solar metallicity; e.g., the simplest case of a 2T 
model with T c = 0.7 keV and Th = 1.5 keV each with the 
same emission measure. Each model has the same total 
emission measure. The spectral energy distribution of the 
IT model has a much narrower peak near 1 keV than the 
2T model; see Figure 5 of Buote (2000a). 

If the observed spectrum is similar to the 2T model, 
such as we observe for NGC 5044 with XMM and Chandra 
within r ~ 30 kpc, then the IT model with solar metal- 
licity is a poor fit to the data. If the iron abundance is 
allowed to vary, then the fitting software (xspec in this 
case) tries to improve x 2 near 1 keV by reducing the value 
of Z-pc while simultaneously raising the continuum to main- 
tain the quality of the fit near the wings of the Fe L region. 
The resulting Fe abundance obtained with the IT model 
is therefore an underestimate of the true value. It is this 
effect that we have previously termed the "Fe Bias" (e.g., 
Buote 2000b,a). 

However, the ability of the \ 2 fitter to increase the con- 
tribution of the continuum within the Fe L region depends 
sensitively on the constraints on the continuum outside the 
Fe L region. At higher energies the model compensates for 
this increase in the continuum by simultaneously decreas- 
ing the abundances of Mg, Si, and S while maintaining an 
approximately fixed a/Ye ratio as discussed in §3 and 4; 
note that Mg, Si, and S have the strongest emission lines 
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Fig. 7. — Radial profiles (units — bottom: arcminutes, top: kpc) of the Fe abundance and associated la errors for the IT (2D) model 
obtained from fitting data with lower energy limit £ m j n = 0.3, 0.5, 0.7 keV. 



at energies above the Fe L region in NGC 5044. 

At energies below the Fe L region (i.e., below 0.7 keV), 
the strongest line is far and away O VIII ly-a near 0.65 keV 
- very close to the bottom edge of the Fe L region. Below 
the O VIII line down to the ~ 0.3 keV bandpass limits of 
the XMM and Chandra CCDs there are no strong emis- 
sion lines. (Allowing the C and N abundances to vary 
affects the inferred value of Z-p e by < 10%.) Therefore, 
possible changes in the continuum level below the O VIII 
ly-a line requested by the \ 2 fitter to better match the 
Fe L region cannot be compensated for by changes in the 
elemental abundances. Changing the temperature is also 
not an option since it is tightly constrained by the shape 
of the continuum both above and below the Fe L region. 
(It is found that the inferred temperature of the IT model 
is quite insensitive to the fitted - or assumed - value of 
Zp e over a large range in Zp e .) Hence, below w 0.7 keV 
the continuum data are a strong inhibitor of the Fe Bias, 
though they do not prevent it entirely as demonstrated by 
the differences we obtained for IT and multitemperature 
models discussed in §2. 

Consequently, it is expected that for spectra where 
the Fe Bias is unimportant (e.g., the XMM spectra for 
r > 40 kpc in NGC 5044) the measured value of Zp e should 
not be very sensitive to the value of E m i n - the lower energy 
limit (0.3-0.7 keV) used in the spectral fitting. In contrast, 
for spectra where the Fe Bias is important (e.g., the XMM 
spectra for r < 30 kpc in NGC 5044) the measured value 
of Zp c should be very sensitive to the value of E m i n , and 
specifically Zp e should decrease as -E m i n increases toward 
0.7 keV. 

This behavior is observed for NGC 5044. In Figure 
7 we plot Z Fc (R) for the IT (2D) model fitted jointly 
to the Chandra and XMM data for different values of 
£ min = 0.3, 0.5, 0.7 keV. For R > 40 kpc the values of Z Fe 
for E m i n — 0.3, 0.7 keV agree within their ~ la errors. The 
values for £" m i n = 0.5 keV are somewhat below these val- 
ues because of the degeneracy with the oxygen abundance. 
(If the oxygen abundance is tied to the iron abundance in 
their solar ratio then the Z Fe values for E m i n = 0.5 keV 



are within 10%-20% of the values of the other -E mm .) For 
R < 30 kpc, where we have the evidence for multitem- 
perature gas and the corresponding Fe Bias, we observe 
that Z Fe decreases as E m i n increases toward 0.7 keV as 
expected. 

We mention that the Zp c values inferred from 2T models 
also follow the same trends, though not as dramatically as 
for the IT models within R < 30 kpc. The reason why 2T 
models are affected at all is that the need for multitemper- 
ature models decreases (though is not removed) as E m i n 
approaches 0.7 keV as discussed in Section 4.4 of Paper 
1. Since a single temperature component dominates more 
in the fits for larger E m { ni the multitemperature models 
suffer more from the Fe Bias for larger E min . For exam- 
ple, the best-fitting values of Z Fe for the 2T (2D) model 
in radial bin #2 are respectively (1.38, 1.19, 0.86) in solar 
units for E min — 0.3, 0.5, 0.7 keV. These values still exceed 
their respective purely IT (2D) counterparts by a factor 
of » 2. 

Hence, for T ~ 1 keV plasmas consisting of multiple 
temperature components distributed across T w 0.7 — 
1.5 keV, we conclude that reliable measurements of Zp c 
can only be obtained for -E nml below w 0.5 keV since the 
continuum emission below 0.6 keV serves as an important 
check on the Fe Bias. In this paper we have decided to em- 
phasize results obtained for E min — 0.5 keV rather than 
for -E m i n = 0.3 keV because of the somewhat better agree- 
ment obtained for the values of Zp e between the XMM and 
Chandra detectors for E min — 0.5 keV discussed in §6.1. 
However, the 10%-20% larger values of Z Fc obtained from 
the joint XMM-Chandra fits for E min = 0.3 keV should be 
considered a reasonable estimate of the systematic error 
arising from the choice of i? m i n for the results quoted in 
this paper. 

6.4. Variable Ah and Intrinsic Absorption 

Allowing for absorption by cold gas (T < 10 4 K) with 
a hydrogen column density (Ah) in excess of the Galactic 
value (A^f 3,1 = 5 x 10 20 cm -2 ) affects the values of Z Fc ob- 
tained from single-temperature models in much the same 
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way as increasing E min toward 0.7 keV with A H = N^ 1 
as discussed above in §6.3. That is, the values of Z Fe ob- 
tained from IT fits with E min = 0.7 kcV and N H = N§ al 
are broadly similar to those obtained from IT fits with 
7V H in excess of TVf al for E min = 0.3 - 0.7 keV; our 
fiducial absorber model is a foreground screen (A(E) = 
exp(—Nna(E))) with separate values of Nn for each annu- 
lus, although we have explored a suite of absorber models 
(e.g., with redshift at NGC 5044) and have obtained fully 
consistent results. 

As discussed in §6.4 of Paper 1 the 2T model is still 
clearly preferred over the IT model when allowing for vari- 
able TVh in each case. Unlike for the IT models, the val- 
ues of Z Fe obtained for the multitemperature models with 
variable Nn only differ significantly with the correspond- 
ing model with Galactic Ah in shell #2; e.g., for the 2T 
(2D) model we have, Z Fc = 0.96 ± 0.05Z Q (variable N H ) 
and Z Fc = 1.19 ± 0.06Z Q (Galactic A H ) - though in both 
cases Z Fc is near the solar value. Since (as discussed in 
Paper 1) the multitemperature models with Ah = A^j* al 
provide better fits than IT models with variable Ah within 
the central ~ 30 kpc, there is no obvious sharp absorption 
feature in the spectrum, and there is no evidence from ob- 
servations in other wavebands for the large quantities of 
cold absorbing material (> 1O 9 M ) implied by the fitted 
values of Ah, we do not take seriously the results obtained 
from the intrinsic cold absorber models. 

We mention that a collisionally ionized "warm" absorber 
(T « 10 5 ~ 6 K) model, in contrast to the cold absorber, 
does not affect significantly energies below ~ 0.5 keV and 
results in fitted values of Z-p c quite similar to the models 
without any intrinsic absorption. Another interesting fea- 
ture of the warm absorber is that the fitted oxygen abun- 
dances in the hot gas in the central regions are w 0.6Z Q 
and the Mg/O ratios are near solar. The near-solar Mg/O 
ratios are in much better agreement with expectations 
from SNe enrichment than the values obtained for all other 
models we explored without a warm absorber. However, 
like the cold absorber there is no evidence for the emission 
in other wavebands implied by the relatively large column 
densities obtained in the X-ray fits for the warm absorber, 
nor is it easy to understand how the temperature of the 
warm gas is maintained. Consequently, we do not discuss 
either the warm or cold intrinsic absorbers further. 

6.5. Background 

We considered the effect of errors in the background nor- 
malization on the measured values of Z Fc . Apart from the 
few emission lines from calibration sources in the XMM 
and Chandra CCDs, the background is a smooth and 
slowly varying function of energy over the 0.3-5 keV band- 
pass. Consequently, if the background contribution to the 
spectrum is underestimated, one will believe the contin- 
uum emission is larger than in reality. As a result, one will 
mistakenly infer smaller equivalent widths for the emission 
lines and therefore underestimate the values for the metal 
abundances. Over-subtraction of the background leads to 
larger equivalent widths and overestimates of the metal 
abundances. 

To examine the sensitivity of the measured values of Z Fe 
to reasonable background errors we repeated our analy- 

4 We note that allowing iVn to vary barely improves the fit: \ 2 



sis using the background templates renormalized to have 
count rates +20% and -20% of their nominal values. By 
fitting IT (2D) and 2T (2D) models to the background- 
subtracted data using these renormalized templates, we 
calculated the variation in Z Fe as a function of radius for 
each case. In the outer annuli (bins #7 and #8), where 
the background is most important, we find that the mea- 
sured values of Z Fc change by < ±15% with respect to the 
nominal background case. The variation is progressively 
smaller at smaller radii being < ±1% in the central bin. 

We mention that even extreme errors in the back- 
ground do not generate qualitatively different results in 
the measured values of Z Fe , especially within the central 
R ~ 30 kpc. For example, if we do not subtract the back- 
ground at all the value of Z Fe is underestimated by < 10% 
for IT (2D) models and overestimated by < 5% for 2T 
(2D) models within radial bins 1-3. The reason why the 
2T models here give a slight overestimated of Z Fc when 
the background is under-subtracted is because the excess 
continuum at higher energies is interpreted as part of the 
higher temperature component. The increased normaliza- 
tion of this component causes an increase in the inferred 
value of Z Fe according to the Fc Bias mechanism (e.g., see 
discussion at the beginning of §6.3). At the largest radius 
(bin 8), if the background is not subtracted at all we ob- 
tain best-fitting values: Z Fc — 0.1 7Z Q compared to the 
nominal value of Z Fe — 0.28Z Q for IT, and Z Fe = 0.33Z Q 
compared to the nominal value of Z Fc = 0.36Z Q for 2T. 

7. EMISSION-WEIGHTED AVERAGE ABUNDANCES AND 
THE FE BIAS REVISITED 

7.1. Central Regions 

Irrespective of whether the hot gas at each radius is 
actually single-phase or multiphase, all models of the spa- 
tially varying spectra that we investigated in Paper 1 re- 
quire that the accumulated spectra within R ~ 30 kpc 
contain a range of temperature components with T « 
0.7 — 1.2 keV. Fitting a IT model to this entire accu- 
mulated spectrum must fail according to the Fe Bias as 
summarized at the beginning of §6.3. 

We demonstrate this manifestation of the Fe Bias as fol- 
lows. We extracted the accumulated spectra of the XMM 
EPIC MOS1, MOS2, and pn data within a radius of 5' 
(48 kpc); we choose a radius slightly larger than 30 kpc 
to fully enclose the peak of the IT temperature profile 
(Paper 1, Figure 3) and to facilitate comparison to pre- 
vious ASCA results below. (Note the Chandra ACIS-S3 
data do not extend to this radius.) The result of fitting a 
IT APEC model simultaneously to the MOS1, MOS2, and 
pn data is shown in Figure 8. Readily apparent are the 
residuals characteristic of the Fe Bias seen in the smaller 
apertures in Paper 1 and in the larger apertures in our 
previous ASCA studies of the brightest elliptical galaxies 
in centrally E-dominated groups (see especially Figure 5 in 
Buote 2000a and the Appendix in Buote 2000b). The fit 
is formally of very poor quality (x 2 = 3354.5 for 816 dof) 
and the iron abundance is 0.51±0.01Z Q . We reiterate that 
the poor fit and low value of Z Fc are exactly as expected 
because of the presence of multitemperature temperature 
components 4 within the aperture - but we have made 

= 2958.1 for 815 dof. The fitted absorption column density is 
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Fig. 8. — MOS1, MOS2, and pn spectra accumulated within a circular aperture of radius, R = 5' (48 kpc), fitted with (Left panel) a single 
temperature (IT) model and (Right panel) a two-temperature (2T) model; both models are 2D. In each case the APEC plasma model is used 
and the solar abundances are taken from Grevesse & Sauval (1998) which use the new (smaller) photosphcric value for the iron abundance. 
The models are the same as discussed in §2 for the spatially resolved analysis. That is, Galactic absorption is assumed and the following 
metal abundances are free parameters: O, Ne, Mg, Si, S, Fe and all other abundances are tied to Fe in their solar ratios. For the 2T model 
the abundances of each temperature component are tied together in the fits. 



Table 3 

Emission- Weighted Average Abundances Within R = 5' = 48 kpc 





Best 


AStatistical 


AModcl 


APlasma Code 


Ai^min 


A7V H 


ABackground 




1.09 


±0.04 


±0.01 


-0.05 


+0.23 


±0.03 


±0.01 


Zsi/ %Fe 


0.832 


±0.017 


±0.004 


+0.058 


+0.013 


±0.011 


±0.003 


Zs/ ZFe 


0.542 


±0.020 


±0.006 


±0.028 


-0.021 


±0.002 


±0.002 


Zol Zp e 


0.331 


±0.028 


±0.005 


-0.052 


+0.089 


±0.014 


±0.003 


Z~Mg/ ZFe 


0.874 


±0.026 


±0.013 


±0.026 


+0.055 


±0.023 


±0.005 


Zn c 


0.56 


±0.09 


±0.10 


-0.37 


+0.26 


±0.04 


±0.03 



Note. — Selected abundances and ratios obtained from fitting the accumulated XMM EPIC MOS1, MOS2, and pn spectral data within a 
circular aperture of radius 5' (48 kpc). The "Best" column indicates the average of the best-fitting values obtained from the 2T and PLDEM 
APEC models expressed in terms of the solar abundance table of Grevesse & Sauval (1998). The "Statistical" column gives the statistical 
error. "Model" is the difference in the best-fitting values obtained from the 2T and PLDEM models. "Plasma Code" represents the difference 
obtained when using the MEKAL plasma code. "-E m i n " is the difference when setting the lower energy limit to 0.3 keV. "iV H " indicates the 
error if the Galactic absorption column density is set to ±10% of the assumed value of 5 X 10 20 cm' 2 . "Background" represents the error if 
the background level is set to ±20% of nominal. Calibration error should contribute less than a 10% error on Zp e considering comparisons to 
ASCA (§7) and Chandra data (§6.1) of NGC 5044. 



no assumption about whether these components 
arise from a radially varying single-phase or a true 
multiphase medium. 

A vastly improved fit which (1) eliminates the residuals 
near 1 keV to a magnitude similar to that present at other 
energies, and (2) provides a value of Zp e over twice as large 
as the IT value is obtained using a simple 2T model (i.e., 
discrete temperature distribution) or a PLDEM (i.e., con- 
tinuous temperature distribution). Each of these models 
adds only 2 free parameters over those of the IT model. In 
Figure 8 we display the 2T model fitted to the total accu- 
mulated EPIC spectra within 48 kpc. It is clear that the 
residuals near 1 keV are greatly reduced the fit is improved 
dramatically (x 2 = 1202.3 for 814 dof). The temperatures 



obtained for the 2T model are T c = 0.804 ± 0.004 keV 
and Th = 1-38 ± 0.02 keV which are very similar to the 
range of values obtained from the spatially resolved IT 
and 2T models (Paper 1); note that the best-fitting ratio 
of emission measures of the hotter and cooler components 
is 1.04. Moreover, the iron abundance for the 2T model is 
a factor of 2.1 times larger than obtained for the IT model: 
Z Fc = 1.08±0.04Z Q (statistical error). The PLDEM gives 
a fit and abundance values extremely similar to the 2T 
model: \ 2 = 1235.1 for 814 dof and Z Fc = 1.10 ± O.O4Z . 
(Note we obtain the following temperature parameters for 
the PLDEM: a = -1.2 ± 0.2, T min = 0.665 ± 0.008 keV, 
T max -T min = 1.12 ± 0.05 keV.) 

The superb agreement between the 2T and PLDEM 



ATVh ~ 5 X 10 20 cm 2 above the Galactic value, implying large amounts of cold gas that have never been observed in other wavebands 
in NGC 5044 or similar systems; i.e., -M a b a Ri m. H A7V H 7r(48kpc) 2 « 3 X 1O 1O M assuming solar abundances and that the absorber is uniformly 
distributed. If the abundance ratios are similar to the hot gas, then Zq/Z Fc 0.2 implies M a t, s ~ lO n M0 because oxygen is the primary 
absorber. If the absorber is not uniformly distributed then the quoted values for M^ s are only lower limits. 
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models demonstrates that the EPIC spectral data accu- 
mulated within 48 kpc have both sufficient sensitivity and 
resolution to unequivocally rule out the IT model (as ex- 
pected) while constraining the DEM well enough so that 
fully consistent values for the abundances are obtained us- 
ing very different multitemperaturc models; i.e., the av- 
erage (emission-weighted) abundances within 50 kpc are 
very well constrained by the EPIC data. 

In Table 3 we list the value of Zp c obtained from the 
multitcmperature fits within 50 kpc and present a detailed 
accounting of the error budget following our discussion in 
§6. It is noteworthy that the largest source of error is E m i n 
which leads to a larger value. The emission-weighted value 
of Zp c is a consistent average of the values within 50 kpc 
obtained from the spatially resolved analysis (Figure 1). 

Also shown in Table 3 are the emission- weighted average 
abundance ratios Zsi/Zp e , Zs/Zp e , Zo/Zp e , and ZM g /Zp e 
obtained within 50 kpc and their associated error budgets. 
These ratios are very tightly constrained and agree with 
the mass-averaged values obtained from the spatially re- 
solved analysis (§3 and 4). We mention that the ZM g /Zp e 
ratio is well constrained, and the only source of error we 
believe could possibly account for its relatively large value 
is incomplete subtraction of the Al calibration lines near 
1.4 keV. However, the small error in the abundance ratios 
associated with reasonable background errors make this 
explanation seem implausible as well. 

The Ne abundance and its error budget are also shown 
in Table 3. Since the key Ne emission line is well hidden 
within the Fe L forest, the value of Z^ e is quite susceptible 
to differences in the temperature model and plasma code. 
In fact, it exhibits by far the largest differences between 
the APEC and MEKAL plasma code. Despite the relatively 
large systematic errors, the value of Z^ e is clearly sub-solar 
with probably values between 0.3-0.5 solar. 

Although we have established that calibration errors 
should contribute at most 20% extra uncertainty in the 
measured value of Zp e (§6.1), we can provide a further cal- 
ibration check by comparing to our previous ASCA studies 
of NGC 5044 (Buote & Fabian 1998; Buote 1999). Buotc 
& Fabian (1998) analyzed the ASCA SIS0 and SIS1 spec- 
tra accumulated within circles of R s=s 5' centered on NGC 
5044. They fitted mekal plasma models over 0.5-5 keV 
where (1) all abundances were tied to iron in their solar 
ratio, (2) the solar abundance table of Anders & Grevesse 
(1989) was used, and (3) the absorption column density 
was a free parameter. They obtained Zp c — 0.25Z Q for 
the IT model (no error given because of the poor fit) and 
Z Fc = 0.62tomZ Q (90% confidence) for the 2T model; 
the 2T model was a clearly superior fit (see Figures 1 
and 5 of Buote & Fabian 1998). If we perform fits to 
the XMM EPIC data in the same region using exactly 
the same models we obtain Zp c = 0.256 ± 0.003Z Q for 
IT and Z Fc = 0.64 ± O.O2Z for 2T (only statistical error 
quoted) in excellent agreement with our previous results 
from ASCA. (Similar agreement is obtained when com- 
paring identical models to Buote 1999.) The excellent 
agreement between XMM and ASCA implies that 
calibration error cannot be a large contributor to 
error in our measurement of Z Fc , certainly < 10%. 

Hence, the XMM data of NGC 5044, whether extracted 
in a large 50 kpc aperture as done in this section or in the 



smaller apertures in previous sections, fully confirm the Fe 
Bias effect, not only for NGC 5044, but by implication also 
for many other bright ellipticals in centrally E-dominated 
groups we examined in previous ASCA studies (Buote & 
Fabian 1998; Buote 1999, 2000b). 

Our results in this section clearly rule out the claim by 
Loewenstein & Mushotzky (2002) and Mushotzky et al. 
(2003) that the hot gas in the central regions of very lu- 
minous galaxies like NGC 5044 must be spatially isother- 
mal. (Our spatially resolved analysis of the temperature 
distribution in Paper 1 indicates that the gas is not single- 
phase as well.) Previous ROSAT and ASCA studies that 
attempted to measure Z-p e using IT models necessarily ob- 
tained values that were biased low (e.g., Awaki et al. 1994; 
Matsushita et al. 1994; Fukazawa et al. 1996; Matsumoto 
et al. 1997; Arimoto et al. 1997; Loewenstein & Mushotzky 
2002; Davis et al. 1999; Finoguenov & Ponman 1999). (We 
note that Matsushita et al. 2000 obtained near-solar values 
for Z-pe in some bright galaxies with ASCA by adding a 
uniform 20% systematic error across the Fe L region. They 
assumed this systematic error arose from calibration error 
or errors in the plasma code rather than from multiple 
temperature components within a large aperture.) 

We mention that although the 2T and PLDEM models 
provide vast improvements in the fits over the IT model, 
their x 2 /dof values are still formally unacceptable. The 
residuals of these fits do not show any obvious features 
and have the same magnitude throughout the bandpass. 
The ratios of model-to-data are within ±5% and appear to 
be largely the result of small inconsistencies between the 
MOS1, MOS2, pn, and ACIS-S detectors. (Formally ac- 
ceptable x 2 values are obtained, e.g., using just the Chan- 
dra data as shown in Table 4 of Paper 1.) However, since 
there are radial variations of the spectral properties within 
the large aperture, we expect that the models fitted to the 
entire large aperture must fail at some level; i.e., the tem- 
perature components of the 2T model shown in Table 1 do 
vary with radius within R = 5'. 

7.2. Outer Regions 

To obtain the average properties at large radius we ex- 
tracted the EPIC spectra within an annulus of R = 5' — 10' 
(48 kpc - 96 kpc) and performed the same analysis just de- 
scribed above for the central region. The results of fitting 
IT and 2T APEC models are shown in Figure 9. At these 
large radii the IT, 2T, and PLDEM models fit quite simi- 
larly: x 2 = 891.9 for 762 dof (IT), X 2 = 840.0 for 760 dof 
(2T), and X 2 = 842.3 for 760 dof (PLDEM). We do not 
see the residuals near 1 keV characteristic of the Fe Bias. 
This is expected since the effective range of temperatures 
within this aperture is only from T « 1.15 — 1.25 keV 
according to the spatially resolved analysis (see in partic- 
ular Figure 3 of Paper 1). This small temperature range 
does not fully span the Fe L region and will not suffer a 
substantial Fe Bias. 

In Table 4 we present the results for the iron abundance 
and abundance ratios and the estimated error budget anal- 
ogously to that for the central region in the previous sec- 
tion. Although the spatially resolved analysis implies a 
small range of temperatures, we do consider the IT results 
in our error budget for the "Model" simply to be conser- 
vative since the \ 2 fits do not obviously distinguish the 



One Temperature 



Two Temperatures 
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Table 4 

Emission-Weighted Average Abundances Within R = 5' — 10' (48 kpc - 96 kpc) 





Best 


AStatistical 


AModel 


APlasma Code 


Ai^min 


A7V H 


ABackground 


Zp e 


0.44 


±0.02 


±0.05 


+0.04 


+0.09 


±0.01 


±0.04 


Zsi/ Zf c 


0.732 


±0.066 


±0.020 


±0.003 


+0.028 


±0.015 


±0.002 




0.55 


±0.09 


±0.01 


+0.01 


-0.02 


±0.003 


±0.12 


Zo/Zpe 


0.96 


±0.11 


±0.02 


-0.11 


+0.18 


±0.04 


±0.10 


Z~Mg/ Zf c 


0.56 


±0.11 


±0.12 


-0.16 


+0.15 


±0.04 


±0.20 


Zn c 


0.20 


±0.11 


±0.18 


-0.11 


+0.24 


±0.04 


±0.17 



Note. — The columns are the same as described in in Table 3 with one exception: "Model" refers to the average value of the IT, 2T, 
and PLDEM APEC models. This average is calculated by first computing the average for the 2T and power-law DEM to get the average 
multitcmperature value. Then this value is averaged with the IT value to obtain the final average value. Hence, "Model" mostly reflects the 
difference between IT and the multitemperature models. 



multitemperature models from the IT models. (We men- 
tion that the 2T model gives T h = 1.33±O.O2Z and T c = 
0.78±0.08Zq with a best-fitting norm^ / norm c — 5.8. The 
PLDEM model gives a = 4.5±0.4, T min = 0.29±0.13 keV, 
T max - T min = 1.12 ±0.11 keV. The DEMs of these models 
are very peaked around a single temperature « 1.35 keV; 
e.g., the multiphase strength (Buote et al. 1999) of the 
best-fitting PLDEM model is ss 0.10 which indicates a 
nearly single-temperature medium.) 

The value of Zp c is tightly constrained to « 40% of the 
average value of Zp c within R = 5'. Notice also that the 
Zsi/Z-F c and Zg/Zp c ratios are consistent with their cor- 
responding values for R = 5' within their 1 — 1.5<r statisti- 
cal errors. The Mg and Ne abundances have rather large 
systematic errors in this outer region, and thus must be 
considered to be consistent with their values for R = 5'. 
However, Zq / Zp c is significantly larger and has a value 
near solar. As shown in §4, the large value of Zq at large 
radius is due entirely to bin #7 which suggests it may not 
be that large. Values of Zq / Zp c at adjacent radii are con- 
sistent with those at smaller radii within the « la errors. 
If the discrepancy is real, it is difficult to understand why 
Zo/Zpe increases with radius while the other ratios are 



consistent with a constant or decreasing profile. Alterna- 
tively, there could be a calibration issue peculiar to the 
measurement of the oxygen abundance. 

8. CONCLUSIONS 

8.1. Iron Abundance and Bias 

One robust conclusion to be drawn from our spatially 
resolved spectral analysis of the XMM data of NGC 5044 
is that the iron abundance drops from a value near solar at 
r « 30 kpc to a value between 0.3-0.4 solar at the largest 
radius probed (« 100 kpc). This radial decrease in Zp e 
is highly significant considering both the statistical and 
systematic errors as illustrated by the emission-weighted 
average abundances obtained from two large apertures 
(§7): Z Fc = 1.09 ± 0.04Z© (statistical) ± 0.05Z© + 
0.18Z Q (systematic) for R = 5' compared to Zp e = 0.44 ± 
0.02Z Q (statistical) ± 0.10Z Q + 0.13Z Q (systematic) for 
R = 5' — 10'. (Note that for ease of presentation we 
have merely added the systematic errors listed in Table 
3. The reader should refer to that table for the explicit 
breakdown of the errors.) There is no evidence that the 
radial decline of Zp e flattens at large radius. The gradient 
in Zp e at large radius measured by XMM is more precise 
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than obtained from previous ROSAT and ASCA studies of 
NGC 5044 (e.g., David et al. 1994; Finoguenov & Ponman 
1999; Buote 2000a). We quote values of Zp e with respect 
to the new solar photospheric value (e.g., McWilliam 1997; 
Grevesse & Sauval 1998) which also agrees with the mete- 
oritic value. 

A second robust conclusion is that the values of Zp c 
within the central regions (r < 30 kpc) are not highly sub- 
solar; i.e., not less than 0.5Zq. In contrast, we find that 
simultaneous spectral fits to the XMM and Chandra data 
within the central regions give values of Zp c ~ IZq for all 
deprojected temperature models assuming only foreground 
Galactic absorption. For the preferred multitemperature 
models (2T and PLDEM) we obtain Z Fc w 1.2Z Q over 
5—30 kpc which dips to Zp c w O.8Z for r < 5 kpc. (These 
values of Zp c obtained from multitemperature models as- 
sume that the abundances are the same for each temper- 
ature component.) The emission- weighted average value 
of Zp e quoted above for R = 5' (48 kpc) rules out with 
high significance a very sub-solar average value of Zp e . 
We conclude that the very sub-solar Fe values obtained 
from previous analyses of ASCA data in the central re- 
gion of NGC 5044 and related systems were affected by 
the Fe Bias as argued previously by us (Buote & Fabian 
1998; Buote 1999, 2000b,a) and others (Allen ct al. 2000; 
Molendi & Gastaldcllo 2001). Recently, we have obtained 
a similar result for the galaxy group NGC 1399 with XMM 
data (Buote 2002). 

We emphasize that for moderate resolution spectra, such 
as provided by the CCDs of the XMM and Chandra satel- 
lites, the Fe abundance of ~ 1 keV systems like NGC 5044 
is not determined from a single isolated emission line, but 
rather from the unresolved Fe L shell emission lines that 
dominate the spectrum between « 0.7 — 1.2 keV. Conse- 
quently, no single Fe emission line is responsible for the 
Fe abundance inferred from current CCD studies. Hence, 
while 20%-30% errors in some Fe L lines still exist in 
the APEC and MEKAL plasma codes (Behar et al. 2001), 
errors of this magnitude cannot be responsible for large 
systematic errors in the inferred Fe abundances from the 
moderate-resolution CCD data. The small (5%-10%) sys- 
tematic offset in the Fe abundances deduced from the 
APEC and MEKAL codes for NGC 5044 (§6.2), each of which 
has different line libraries and approximation methods, in- 
dicates the magnitude of the systematic error in the Fe 
abundance expected from remaining errors in the model- 
ing of the Fe L lines. 

Of more serious concern is that Fe abundance measure- 
ments can be systematically biased by improper definition 
of the continuum, inclusion of an intrinsic absorber com- 
ponent when none exists, and forcing a single-temperature 
component to fit a multitemperature spectrum possessing 
temperature components near 1 keV. So long as the lower 
energy range used for spectral fitting is well below the Fe 
L region (i.e., at least as low as 0.5 keV), the bias arising 
from a poor continuum definition is mitigated (§6.3). Since 
IT models with intrinsic absorption do not fit as well as 
2T and PLDEM models with Galactic Nn, and there is no 
evidence outside the X-ray band for large quantities of cold 
(or warm) absorbing gas in NGC 5044 or other X-ray lu- 
minous groups and clusters, it is reasonable to include only 
foreground absorption from the Milky Way (§6.4). Finally, 



the Fe bias associated with forcing a single-temperature 
component to fit the multitemperature spectra in the cen- 
tral regions of NGC 5044 results in underestimates of Zp e 
often by as much as 50% or more. 

8.2. Implications for Supernova Enrichment 
8.2.1. Iron 

The approximately solar values for Zp e within the cen- 
tral ~ 30 kpc determined from our X-ray observations of 
the hot gas exceed the typical values of « 0.5Z Q obtained 
for the stars averaged over a region of half an effective ra- 
dius (Trager et al. 2000). This result implies that the hot 
gas in NGC 5044 could have been enriched by a substantial 
amount of Type la supernovae (SNe la) which alleviates 
the most serious conflict of previous X-ray observations of 
elliptical galaxies and centrally E-dominated groups with 
the models of chemical enrichment that assume a Milky 
Way IMF and SNe la rate as observed in elliptical galax- 
ies (Cappellaro & Turatto 2002). 

This conflict, called the "Iron Discrepancy" in elliptical 
galaxies by Arimoto et al. (1997) and Renzini (1997), is 
partially resolved for NGC 5044 by our XMM and Chan- 
dra observations which demonstrate that previous under- 
estimates of the iron abundances with ASCA arose par- 
tially from the Fe Bias and partially from the use of the 
incorrect solar abundance standard. Our recent XMM ob- 
servation of the centrally E-dominated group NGC 1399 
(Buote 2002) provides similar results for the iron abun- 
dance. Other XMM observations (EPIC and RGS) of 
dominant ellipticals in groups (e.g., NGC 533, Peterson & 
et. al. 2003; NGC 4636, Xu et al. 2002) and clusters (e.g., 
M87, Molendi & Gastaldello 2001; Gastaldello & Molendi 
2002; A2029, Lewis et al. 2002) find near solar central Fe 
abundances consistent with significant SNe la enrichment. 

8.2.2. a/Fe Ratios 

The ratios of the abundances of the elements O, Mg, 
Si, and S to Fe provide strong constraints on the frac- 
tion of iron in the hot gas produced by SNe la. As we 
discussed in §3 and §7, next to Fe the best constrained 
abundance is that of Si. The emission-weighted aver- 
age Zsi/ Zp e ratio within R = 5' (48 kpc) is, Zs;/^fo = 
0.832 ± 0.017 (statistical) ± 0.018 + 0.071 (systematic) in 
solar units. This ratio translates to a SNe la fraction of 
67%- 79% where we have used (1) the la statistical error 
on Zsi/Zp e , (2) the theoretical SNe la yield of the "con- 
vective deflagration" W7 model of Nomoto et al. (1997b), 
and (3) the range of theoretical SNe II yields reported in 
Gibson ct al. (1997). Note that for this calculation we 
have used the abundance ratios for the theoretical super- 
nova models expressed in terms of the Grevesse & Sauval 
(1998) solar abundance standard as reported in Table 2 
of Gastaldello & Molendi (2002). Using the la range for 
Zsi/Zp e at large radius (i.e., for R = 46 — 96 kpc) we 
obtain a SNe la fraction of 70%-86%, consistent with the 
result within R — 48 kpc. 

The emission- weighted average sulfur-to-iron abundance 
ratio within R = 5' (48 kpc) is, Z s /Z Fc = 0.542 ± 
0.020 (statistical) ± 0.010 + 0.007 (systematic) in solar 
units, which is consistent with the average value obtained 
at larger radii out to 100 kpc. Using the la statistical er- 
ror on the X-ray measurement, along with the theoretical 



15 



SNe models mentioned above, we obtain a SNe la fraction 
of 52%-82%, in good agreement with that inferred from 
the Zsi/Zf c ratio. 

The largest contributor to the uncertainty in the SNe la 
fraction deduced from these abundance ratios is the theo- 
retical SNe yields. For example, if we assume the SNe II 
yields of Nomoto et al. (1997a) we obtain SNe la fractions 
of 76%-77% from Z Si /Z Fc and 70%-75% from Z s /Z Fc us- 
ing the la statistical errors on the X-ray measurements 
within R = 48 kpc. The so-called "delayed-detonation" 
SNe la models, WDD1 and WDD2, of' Nomoto et al. 
(1997b) appear to be incompatible with the measured 
Zgi/Zpe and Zs/Z Fe ratios. The WDD1 model predicts 
Zsi/ Z Fo = 1.15 and Zq/ Z Fc = 1.01 in solar units which are 
incompatible with our measurements. The WDD2 model 
predicts Zs/Z Fc = 0.62 which is Aa above the measured 
value within R = 48 kpc. However, if the various sys- 
tematic errors shown in Table 3 are considered then the 
discrepancy is not highly significant. Nevertheless, the 
Zsi/Z Fe , and especially Zg/Z Fe , ratios appear to prefer the 
convective deflagration W7 SNe la model over the delayed- 
detonation models. 

Hence, we infer that SNe la have contributed 70% — 80% 
of the iron mass within a 100 kpc radius of NGC 5044 - 
a result that is fully consistent with that obtained from 
our recent analysis of an XMM observation of the bright 
nearby galaxy group, NGC 1399 (Buote 2002). This 
SNe la fraction is also similar to that inferred for the Sun 
and therefore suggests a stellar initial mass function sim- 
ilar to the Milky Way (e.g., Rcnzini et al. 1993; Rcnzini 
1997; Wyse 1997). 

We mention that the uncertain Ne abundance is also 
consistent with this result - see §7. The average Zo/Z Fe 
ratio (except bin #7) is also consistent with this level of 
SNe la enrichment: 75%-89% using the la statistical er- 
ror on the X-ray measurement within R — 48 kpc, the W7 
SNe la model, and the range in SNe II yields from Gib- 
son et al. (1997). Like Zo/Z Fc in bin #7, the near-solar 
Mg/Fe ratio in the central regions implies a much smaller 
SNe la fraction (ss 50%). We cannot explain the apparent 
anomalous behavior of these Zo/Z Fe and ZMg/Z Fe ratios, 
but since their relevant lines are blended with strong Fe L 
lines and (in the case of Mg) calibration lines, we consider 
these discrepancies to be tentative. 

8.3. Central Iron Deficit and Problems with Enrichment 

Models 

Although the near-solar central value of the iron abun- 
dance and the values of the a/Fe ratios provide strong ev- 
idence for SNe la enrichment of the hot gas from the cen- 



tral galaxy, there remain some problems with this scenario. 
First, the abundances in the hot gas generally rise toward 
the center and then dip in the innermost bin; this dip 
is more pronounced for the deprojected abundance pro- 
files. Such a dip is difficult to understand in the context of 
enrichment from the central galaxy which should be most 
pronounced at the very center. Some cD clusters also show 
such a dip (e.g., Centaurus, Sanders & Fabian 2002) which 
Morris & Fabian (2002) have suggested is an artifact of 
attempting to model a highly inhomogeneous metal distri- 
bution with a homogeneous spectral model. This model is 
attractive since it could explain why sub-keV cooling gas 
has not been found by Chandra or XMM in cooling flows. 

In NGC 5044 we find that the central abundance dip 
is indeed sensitive to assumptions in the spectral model. 
(Central dips are also sensitive to the lower energy limit 
of the bandpass as discussed in §6.3.) If we tic oxygen 
to iron in their solar ratio the dip in the abundances is 
less pronounced. In the spirit of Morris & Fabian (2002), 
if instead we allow the abundances on each temperature 
component of the 2T model to vary separately then the 
abundances on the cool component remain large all the 
way to the center. We show the central iron abundance 
profile of this model in Figure 10. Because of the large 
number of free parameters in this model the precise ra- 
dial variation of Z Fc in Figure 10 should be treated with 
caution. But the basic idea that the property of a central 
dip is sensitive to further assumptions about the spectral 
model is clear. 

Second, we note that although correcting for the Fe Bias 
partially removes the "Iron Discrepancy" noted by Ari- 
moto et al. (1997), chemical models of elliptical galaxies 
without cooling flows predict central iron abundances even 
larger than we have measured for N5044 and NGC 1399 
(Brighenti & Mathews 2003). 

It should also be emphasized that the properties of the 
central ~ 5 kpc (i.e., central R s=s 30") are distinctly differ- 
ent from those at immediately larger radii. Hence, studies 
of NGC 5044 with the high-resolution gratings of Chandra 
and XMM can only probe this special region and cannot 
tell us about the properties at larger radii. We note, how- 
ever, that the RGS, EPIC, and ACIS-S3 data give consis- 
tent results in their overlap regions (§6.1). 

It is a pleasure to thank T. Fang for communicating 
results of his unpublished analysis of the RGS data of 
NGC 5044. We gratefully acknowledge partial support 
from NASA grants NAG5-9956, NAG5-10758, and NAG5- 
10748. 
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